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Abstract. We present photometric and spectroscopic observations of the late afterglow of GRB 000911, starting 
~ 1 days after the burst event and lasting ~ 8 weeks. We detect a moderately significant re-brightening in the 
R, I and J lightcurves, associated with a sizable reddening of the spectrum. This can be explained through 
the presence of an underlying supernova, outshining the afterglow ~ 30 days after the burst event. Alternative 
explanations are discussed. 
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1. Introduction 

An anomalous re-brightening has been detected in the 
optical lightcurves of at least two gamma-ray bursts 
~ 30 days after the burst event (GRB 980326: Bloom 
et al. 1999; GRB 970228: Reichart 1999, Galama et al. 
2000). These have been tentatively interpreted as due to 
the simultaneous explosion of a supernova (SN) with a 
lightcurve similar to that of SN1998bw (Galama et al. 
1998) which outshines the afterglow emission at the time 
of the SN peak. GRB 980326 showed a clear excess in 
R band (aided by having a very faint host galaxy) and 
the spectrum changed from blue (early) to red close to 
the peak of rebrightening; however, neither the redshift of 
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the afterglow nor the host is known to date (Bloom et al. 
1999). GRB 970228 had the advantage of a known redshift 
(Djorgovskij et al. 1999), a moderate amount of multiband 
data but lacked spectroscopic coverge. A broad band mul- 
tifilter spectral energy distribution (SED) of GRB 970228 
at day ~ 30 was extracted (Reichart 1999), but the K 
band measurement, important to constrain the reddening, 
was questioned by Galama et al. (2000). 

The confirmation of a GRB-SN connection is clearly 
of primary importance to solve the problem of the burst 
progenitor (Eichler et al. 1989; Woosley 1993; Paczyhski 
1998), but also bears important consequences on the 7-ray 
radiation mechanism (Lazzati et al. 2000) and the origin 
of emission and absorption features in the X-ray spectra 
of the burst proper and its afterglow (Lazzati et al. 1999; 
Amati et al. 2000; Piro et al. 2000; Rees & Meszaros 2000; 
Vietri et al. 2001; Lazzati et al. 2001). 

Here we present the results of a simultaneous multi- 
filter observational campaign designed at detecting and 
studying the spectrum of the re-brightening component 
in the burst of September 11 th , 2000. We observed the 
optical transient (OT) in five filters at three epochs. A 
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low resolution spectrum was taken ~ 36 days after the 
burst explosion. 

GRB 000911 was detected by Ulysses, NEAR and 
Konus on September 11.30237 (Hurley 2000). It had a 
duration of ~ 500 s, a 25-100 keV fluence of T ~ 
5 x 10~ 6 erg cm' 2 and a peak flux over 0.5 s of ~ 
9 x 10~ 7 erg cm -2 s _1 . Its afterglow was detected at 
radio wavelengths (8.46 GHz, Berger & Frail 2000) and 
subsequently confirmed in the optical R band (Berger et 
al. 2000) at a position of R.A. 02 h 18 m 34 s .36 and Dec. 
07°44'27".7 (Price 2000). The OT displayed a fairly rapid 
decay law F v [t) oc t~ (Price et al. 2000, see also 



Sect. 3.1 ) in the R band. The redshift of the host galaxy 
was measured to be z = 1.06 (Price et al. 2001). 

2. Observations 

We observed the OT associated with GRB 000911 with 
the ESO/VLT-Antu telescope (instruments, FORS1 and 
ISAAC), the Keck I telescope (instrument, LRIS) and the 
MSO 50-inch telescope; see Tables 1 and 2 for log of ob- 
servations. The VLT observations were designed to opti- 
mize the chances of detecting a re-brightening due to a SN 
similar to that of SN 1998bw. Furthermore, we obtained 
low resolution spectra with the FORS1 instrument (grism 
1501, blocking filter OG590) on day 36, around the time 
of the expected peak SN emission. 

2.1. Photometry 

All scientific frames were reduced in a standard way with 
the ESO-MIDAS package (98NOV and 99NOV versions). 
Photometric analysis was carried out with DAOPHOTII 
(Stetson 1987) as implemented in the ESO-MIDAS and 
IRAF (version 2.11.3) packages. We derived both aper- 
ture and point-spread function (PSF) fitting photometry. 
The model PSF was always based on a large number of 
stars. Seeing conditions were good in the first and third 
VLT runs, (~ 0.8arcsec), and exceptional in the second 
one (~ 0.5arcsec). For all three VLT runs suitable stan- 
dard star fields were observed: SA98, Mark A and SA113 
in the optical; FS28, S234 and P530 in the J band. The 
observations were all performed under photometric con- 
ditions. A comparison of our photometric calibration and 
that by Henden (2000) shows that ours is slightly fainter 
(~ 0.1 mag) for stars at the faint end of Henden's list. 

The observation log and photometric measurements 
for the three multifilter VLT observations are reported 
in Tab. [I]. Table shows the additional photometric mea- 
surements taken with the Keck and MSO50 telescopes, as 
well as two photometric points extracted from the VLT 
spectrum. All these measurements, together with two ad- 
ditional photometric points taken from the literature (see 
Tab. ||), are shown in Fig. 0. The OT appears point like in 
all filters at all times, with the exception of the B and V 
images at day 53, in which the OT is not consistent with 
the PSF of bright stars. The faintness of the OT in B and 
V bands at day 53, however, makes it impossible to mea- 



sure a reliable extension and ellipticity of the source. The 
error bars of photometric points were derived following the 
standard procedure. All pertinent factors affecting photo- 
metric precision are taken into account: statistical errors 
on the measured counts converted to electron units, CCD 
read-out noise, calibration errors applying the usual error 
propagation formulas. On the other hand, the photome- 
try of the OT was not affected by the galaxy background 
possibly present in our images. In fact, we did not detect 
significant evidence of deviation from the field star PSF for 
the OT (the only exception is for the last B and V images, 
see above) . The underlying galaxy within the limit of our 
observational data seem to be consistent wih a point -like 
shape perfectly superposed (i.e. within a fraction of the 
FWHM of our images) to the OT. After having removed 
all star-like objects in our images no diffused background 
emerged in the OT region. A systematic uncertainty of 
0.05 mag was added in quadrature to all data, in order 
to get rid of the differences possibly introduced by images 
taken with different telescopes. 

2.2. Spectroscopy 

The spectral frame reduction was again performed in a 
standard way with ESO-MIDAS tools (99NOV version). 
The extraction of the scientific spectra was then obtained 
with IDL routines written by us. Observing conditions 
were photometric. However the closeness of the moon to 
the target (~ 30°), the seeing variations (from ~ 0.8 arcsec 
to ~ 1.8 arcsec) and the OT faintness made the spectrum 
extraction critical. In particular it proved necessary to re- 
bin our spectra to the much lower resolution of 110 A per 
pixel in order to increase the S/N. Two independent spec- 
tra were extracted from different frames, one hour expo- 
sure each, and found to be consistent with each other. 
Absolute calibrations were performed through the obser- 
vation of a standard star (EG21) with the same instru- 
mental set-up. 

The resulting average spectrum is shown by the filled 
circles in Fig. ||. The spectrum, despite its very low res- 
olution, was determined to be considerably redder than 
any afterglow spectrum observed so far. A power-law fit 
yielded F(v) oc ^- 5 - 3±0 - 8 (\aerror, \ 2 — 22 for 21 degrees 
of freedom, hereafter d.o.f.). 

3. Data modelling 

We have modeled the data with a composite spectrum 
given by combining an external shock synchrotron com- 
ponent (Meszaros & Rees 1997) plus a host galaxy. In 
addition, we examined the possible role of a supernova 
component. The details of the photometry and spectral 
fitting are discussed in the following subsections. 

3.1. Photometry 

The photometric data were dereddened for Galactic ex- 
tinction according to the maps of Dickey and Lockman 
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Fig. 1. Lightcurves of the afterglow of GRB 000911. From top left to bottom right, the J, /, R, V and B lightcurves are 
plotted. The thick solid curves show the best fits obtained with our three component model. The dashed, dotted and 
dot-dashed lines show the SN, galaxy and ES components, respectively. The thin solid lines in the J, I and R panels 
indicate the best fit for a model comprising only the galaxy and ES (without SN). The thin line is indistinguishable 
from the thick solid line for the V and B filters. 



(1990) (Eb-v = 0.112) and following extinction curves 
from Cardelli ct al. (1989), and converted to flux densi- 
ties. Extinction in the host galaxy was not modeled since 
no additional extinction was required by the data. 

The lightcurve modelling was obtained as the sum of 
the above three contributions. First we considered an ex- 
ternal shock afterglow component (hereafter ES) of the 
form: 

F(u,t) =A ES v~ a t~ 2a . (1) 

This equation holds for a jet geometry after the break 
time (Sari, Piran & Halpern 1999). Such a configuration 
is obtained by a broadband fitting of the GRB 000911 
afterglow (Price et al. 2001). The constant flux from a host 
galaxy was added in the five bands as a free parameter 
(Gb, Gv, Gr, Gi and Gj). 

First, an ES plus galaxy model was fitted to the data. 
The best fit gave a decay slope a = 0.724±0.006 (temporal 
slope 5 = 1.45±0.012) with X 2 = 24.4 for 18 d.o.f.. This is 
an acceptable fit, with chance probability of P ~ 10% to 
obtain a higher \ 2 value. However, the fit can be improved 
by adding the SN component at the redshift z — 1.06 of 



the host galaxy (see Sect. [I]). The lightcurve of the su- 
pernova component was obtained by spline interpolation 
of the data of SN1998bw (z=0.0085; Galama et al. 1998). 
Cosmological parameters Hq = 65 km s _1 Mpc -1 and 
qo = 0.5 were adopted to compute the flux as a function 
of redshift, and the time profile was stretched by a fac- 
tor 1 + z. Following Bloom et al. (1999), the spectrum 
of the supernova was analytically extended in the rest 
frame ultraviolet assuming a power-law F{v) oc v~ 2 - 8 for 
A < 3600 A. SN1998bw was adopted as a template since 
it is the SN with the most widely accepted association 
with a GRB (GRB 980425) due to the coincident position 
(Galama et al. 1998) and its peculiar relativistic expan- 
sion (Kulkarni et al. 1998). In addition, it has a very well 
sampled UBVRI lightcurve (Galama et al. 1998). 

The addition of the supernova component changes 
slightly the spectral slope of the ES component (a = 
0.748 ± 0.006) with X 2 = 13.9 (P - 67%). In principle, 
the addition of this component does not introduce any 
new free parameter since the redshift is known and the 
SN luminosity fixed to the value of SN1998bw. However, 
the lightcurves of different SNe, even within the same 
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Table 1. Log of the VLT photometric observations of the afterglow of GRB 000911. 



At (d) 


Filter 


Telescope 


Exposure time (s) 


Magnitude 


Seeing 


1.435 


R 


MSO50 


1200 


20.70 ± 0.08 


2.1" 


1.435 


V 


MSO50 
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20.90 ±0.1 


2.1" 
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Keck 
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22.30 ± 0.06 


0.6" 
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21.82 ± 0.06 


0.7" 
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14.15 
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VLT Spectrum 
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Table 2. Additional photometric points from Keck, MSO 50 inch telescope and VLT spectroscopy. 



At (d) 


Filter 


Telescope 


Magnitude 


Reference 


Seeing 


111.9 


R 


Keck 


25.4 ±0.1 


Price et al. 2001 


0.7" 


112.0 


I 


Keck 


24.4 ±0.1 


Price et al. 2001 


0.65" 



Table 3. Photometric points from literature. 



class, are different, with different peak luminosity and 
peak time. We hence let the normalization of the SN com- 
ponent free to vary, with the redshift, the color and the 
time profile held fixed. We obtain a normalization factor 
A = 0.9 ± 0.3, with X 2 = 13-8 for 17 d.o.f.. An F test 
applied gives a statistical confidence of 99.8% (2.9a) for 
the fit improvement. Alternatively, we allowed the red- 
shift to vary. Interestingly, we obtain zsn = bl, m good 
agreement with the z — 1.06 measured spectroscopically 
by Price et al. (2001) for the host galaxy of GRB 000911. 

The re-brightening component in the lightcurve is 
hence remarkably similar in luminosity, shape and color 
to the lightcurve of SN1998bw. As a final test, we allowed 
for a temporal shift At between the SN and GRB explo- 
sions. Keeping all the other SN parameter fixed to the 
values of SN1998bw and the redshift z = 1.06, we obtain 
At = 0^7 5 days (1 a), showing that the SN explosion may 
anticipate the GRB but only by ~ 1 week (see, e.g., the 
Supranova model by Vietri & Stella 1998). 

The best fit three component model (ES plus galaxy 
plus SN) is shown by the thick solid line in Fig. [y. The 
figure also shows the three individual components (ES, 
galaxy and SN). 



In Fig. g we report the evolution of the color index 
R — I as a function of time. The color significantly evolves 
from blue to red, as time goes on. A formal best fit with 
a constant color yields \ 2 — 23 for 5 d.o.f.. The reddening 
behavior can be due either to the presence of an underlying 
supernova or to a red host galaxy (R — I ~ 1) which 
certainly dominates the color evolution at late time (t ~ 
100 d). Again, the quality of the data, especially at late 
times, does not allow us to unambiguously select a model. 

3.2. Spectroscopy 

According to the fits described above, none of the three 
components clearly dominated the emission at the time 
of the spectroscopic observation (~ 36 days; see Fig. |l]). 
Moreover, the wavelength range over which the spectrum 
has a sufficient signal to noise to be safely extracted cov- 
ers only a limited spectral range (6000A < A < 8500A). 
For this reason, the spectrum by itself can be modeled in 
different ways, all consistent with the data. 

To better constrain the models, we added photometric 
information to the spectrum. In the upper panel of Fig. [?[ 
the spectrum (filled dots) is plotted together with the J 
band measurement at day 25 (diamond). Note that the 
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Fig. 2. Evolution of the R — I color index with the time 
elapsed after the GRB event. The color shows a significant 
evolution from the value R — I ~ 0.55 in the first 20 days 
to the final value of R — I ~ 0.85 afterwards. 



change in the J magnitude from day 25 to day 36 is ex- 
pected to be small (see the first panel of Fig. |l|). In the 
lower panel, the best fit galaxy components obtained with 
the lightcurve fitting procedure (from left to right Gb, 
Gy, Gr, Gj and Gj) are shown for the five photometric 
filters used (triangles). We modeled the ensemble of these 
data with a galaxy spectrum (templates from Calzetti et 
al. 1994) plus a power-law ES spectrum (parameters fixed 
to the best-fit values from the lightcurve) and a type Ic 
supernova spectrum a few days after the peak (SN 1987M, 
Filippenko et al. 1990)0. 

The two models were fitted as follows: the total tem- 
plate spectrum was fitted to the spectral data together 
with the J band measurement; galaxy magnitudes derived 
from the template galaxy were fitted to the data in the 
lower panel. For the galaxy plus ES model, a formal best 
fit was obtained with a dust-enshrouded starburst galaxy 
template. The fit gave X 2 = 58 for 27 d.o.f. (P ~ 0.05%). 
A better fit { X 2 = 37 for 27 d.o.f., P ~ 9%) was ob- 
tained by adding a SN component, with a moderately 
dust-enshrouded starburst galaxy (0.11 < Eb-v < 0.21) 
as a template. The fit is shown in Fig. || overlaid on the 
data. 

In order to obtain a single statistical indicator combin- 
ing the photometric and spectral information, we finally 
fitted simultaneously all the available data with the appro- 
priate galaxy template. This yielded \ 2 — 44 (39 d.o.f.) 



and 



X 



69 (40 d.o.f.) for the models with and with- 



out SN component, respectively. The x 2 decrease has a 
statistical significance of 4a, according to the F-test. 



1 A spectrum of SN1998bw would be more suited for this 
comparison, however no spectra of SN1998bw at peak extend 
to wavelengths A < 3500 A, a range that is necessary for a 
proper comparison with the data of a SN at Z = 1.06. 



4. Alternatives 

We detected a possible re-brightening and reddening of 
the OT associated with GRB 000911 ~ 30 days after the 
GRB event. We showed that this is consistent with the 
presence of an underlying SN component. In the following 
we discuss two possible alternatives. 

4.1. Dust echoes 

It has been suggested that the re-brightening detected 
in the late optical lightcurves of several GRBs may 
be explained by dust scattering or reprocessing of the 
burst and/or afterglow photons (Waxman & Draine 2000; 
Esin & Blandford 2000). The spectrum and lightcurve 
of these components have been modeled in detail by 
Reichart (2001). He concludes that the re-brightening in 
the lightcurve of GRB 970228 cannot be explained with 
dust scattering and reprocessing, because it predicts a 
spectral shape different from the observed one. 

The same conclusion can be drawn in the case of 
GRB 000911. In fact, the re-brightening due to dust 
reprocessing has an exponential cut-off at wavelengths 
A < 2.2(1 + z) fxm (Reichart 2001), which lies in the 
Johnson L filter for a z = 1.06 burst. We detect a re- 
brightening in the R, I and J filter, and hence this in- 
terpretation can be ruled out. Scattering from dust, on 
the other hand, does not change the spectral index of the 
afterglow by more than ^0.5 for intermcdiatc-low dust 
opacities (Reichart 2001). The spectrum we measure at 
day 36 implies a much larger reddening (see Sect. 2.2) and 
a dust echo for a low-intermediate dust opacity Ta U st ~ 12 
cannot be fitted to the data even considering the contri- 
bution from a red host galaxy (the fit yields a x 2 compa- 
rable to that obtained for the ES plus galaxy model). For 
a Tdust > 12 dust cloud, the re-brightening would be so 
small that it would remain undetected. The above conclu- 
sion can be clearly further strengthened through the mea- 
surement of the J magnitude of the host galaxy, which 
would allow us to better constrain the temporal and spec- 
tral shape of the re-brightcning component. 

4.2. Wind-Fireball Interaction 

In the collapsar model for GRBs (MacFadyen & Woosley 
1999), or in any other model involving a massive star, 
the key to obtain relativistic motions is the escape of an 
energy-loaded fireball from the stellar environment. A fo- 
cused low-entropy jet that has broken free of its stellar 
cocoon is likely to arise from a Wolf-Rayet (WR) progen- 
itor. WR stars are characterized by strong stellar winds 
gradually shedding most of the star envelope. The decel- 
eration of this pre-collapse wind by the pressure of the 
surrounding medium causes the ejected mass to accumu- 
late at a radius ~ 10 17 M_4n^ 2 v 3 *' 2 cm (where M is the 
mass loss rate in units of solar masses per year, rig is the 
density of the surrounding medium in units of cm -3 , v 
the wind velocity in units of km s -1 and we adopt the 
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Fig. 3. Spectrum of the OT of GRB 000911 observed 36 
days after the burst explosion (upper panel, filled dots). 
The spectrum is modelled with a SN type Ic spectrum 
plus a background starburst galaxy and an ES compo- 
nent (see text). The solid line shows the total spectrum 
smoothed with a 110 A boxcar filter, while the dashed 
line represents the supernova component (SN1987M, see 
text). The lower panel shows the template spectrum of 
the best fit galaxy model (thin solid line). Triangles are 
the galaxy photometric measurements as derived from the 
multiband fitting (see sect. 



3.1 



and the dotted horizontal 
lines in Fig. ^). The vertical position of the horizontal bars 
indicate the BVRIJ filter fluxes derived from the galaxy 
template; their width is equal to the full width at half 
maximum of the filters. 

convention Q = 10 x Q x ). The impact between the for- 
ward shock and these high-density regions should be ob- 
served as a re-brightening of the afterglow, typically with 
a redder spectrum. Using detailed stellar tracks for the 
evolution of massive stars, as described in Ramirez-Ruiz 
et al. (2001), we found that the re-brightening of GRB 
000911 at z = 1.06 can be explained by the interaction 
of a relativistic blast wave expanding into the ambient 
medium expected at the end of the life of a 40 Mq (7 Mq 
core) WR star evolving with solar metallicity. The shock 
front expands within a n(r) = 1.5 x 10 35 r~ 2 cm" 1 stellar 
wind until it reaches the density enhancement (see Fig. 3 
of Ramirez-Ruiz et al. 2001). If we consider the forward 
shock emission alone (shock wave into the accumulated 
wind), it is not possible to have a spectral slope steeper 
than a ~ 3.5 (Ramirez-Ruiz et al. 2001). Including the 
emission of the reverse shock, spectral slopes as steep as 
a ~ 4 can be obtained, depending on the radiative effi- 
ciency of the reverse shock. The spectroscopic observation 
can then be reproduced if some of the reddening is at- 
tributed to the host galaxy emission. 

We have fitted this model to the photometric data only 
obtaining a \ 2 intermediate between the fit with a sim- 



ple ES plus galaxy model and the fit with the SN addi- 
tional component. This model also requires the addition 
of 4 free parameters (the Lorentz factor of the fireball, 
the mass of the progenitor star, its metallicity and the 
density of the surrounding medium) and the need for an 
enhanced reverse shock emission to reproduced the ob- 
served spectrum. Since the addition of more degrees of 
freedom does not improve the statistics of the \ 2 1 we con- 
sider the SN a better explanation for the re-brightening 
in the GRB 000911 lightcurve. 

5. Conclusions 

We presented late time multifilter observations of the op- 
tical transient associated to GRB 000911. This set of ob- 
servations was designed to detect and analyze the re- 
brightening associated with (some) GRB afterglows ap- 
proximately one month after their explosion (Bloom et al. 
1999; Reichart 1999; Galama et al. 2000). In addition to 
photometric data, a low resolution spectrum was taken 
~ 36 days after the burst explosion. 

The lightcurve and spectrum were fitted with an ex- 
ternal shock plus galaxy model, with the possible addition 
of a supernova, similar to SN1998bw. The addition of the 
SN component gives a better fit, with a statistical signifi- 
cance of 4 a. Even though the observations were designed 
to detect unambiguously the additional supernova com- 
ponent, the signal to noise ratio that could be achieved 
was limited, due to the combination of two effects. Firstly 
the burst was at a redshift larger than in the two previous 
cases (GRB 970228: z = 0.695; GRB 980326: z ~ 0.9); sec- 
ondly, the host galaxy associated with GRB 000911 had 
a brightness comparable to that of the SN component. 
For this reason, a fourth multifilter observation at a much 
longer time (> 100 days) is mandatory in order to disen- 
tangle the two components. We were able to find R and 
I magnitudes at this late times in the literature (Price et 
al. 2001), but the most important measurement in the J 
band (see the difference between the thick and thin solid 
lines in the J panel of Fig. was not available. 

With the present data, a word of caution should be 
spent, since is not possible to unambiguously assess the 
presence of the SN. However, if future observations will 
allow us to better constraint the magnitude of the host 
galaxy and will confirm the presence of the rebrightening, 
we will be able to disentangle the SN component and to 
provide simultaneous multiband SEDs at the time of our 
3 VLT observations. Such time resolved broad band SEDs 
will allow to better understand the spectral evolution of 
the bump in the lightcurve and hence to understand its 
physical origin. 
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